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Abstract 

We consider a two reservoir model of quantum error correction with a hot bath causing errors 
in the qubits and a cold bath cooling the ancilla qubits to a fiducial state. We consider error 
correction protocols both with and without measurement of the ancilla state. The error correction 
acts as a kind of refrigeration process to maintain the data qubits in a low entropy state by 
periodically moving the entropy to the ancilla qubits and then to the cold reservoir. We quantify 
the performance of the error correction as a function of the reservoir temperatures and cooling rate 
by means of the fidelity and the residual entropy of the data qubits. We also make a comparison 
with the continuous quantum error correction model of Sarovar and Milburn [Phys. Rev. A 72 
012306]. 

PACS numbers: 03.67.Pp, 03.65.Yz 
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I. INTRODUCTION 



It is clear by now that the operation of a large scale quantum computer will require 
some form of protection of the data qubits from decoherence arising from interaction with 
their environment and imperfect control fields which are applied to carry out quantum 
logic gates. There are three methods known to provide such protection — quantum error 
correcting codes DQDDDDDD , decoherence free subspaces |9|, llOL I 111 . Il2l . Il3l | , and 
dynamical suppression of decoherence [141 ] . 

In this paper we focus on the role of ancilla preparation in a quantum error correction 
protocol. We consider a situation where the ancilla qubits are reset to a fiducial state by 
coupling them to a cold reservoir. We model the source of errors as a hot reservoir coupling 
to both the data and ancilla qubits. This system can be viewed as a kind of refrigerator 
where the data qubits are preserved in a low entropy state via the error correction protocol 
with the excess entropy being dumped into the cold reservoir when the ancilla qubits are 
reset. The performance of the code then depends critically upon the rate at which this cycle 
can be performed. The resetting of the ancilla qubits also requires an energy input of at 
least kgT log2 per bit of information erased where T is the temperature of the environment 
into which the information is lost jisl . ^| . This defines a minimum energy cost for running 
the error correction protocol which is different from that considered by one of the authors 



in Ref. [17 1. 

In a real system it may not be practical to control the coupling of the ancilla qubits to a 
cold reservoir at will. It may be necessary to use an active cooling method which could be 
similar to that analyzed in [18[. Nevertheless we model the cooling of the ancilla qubits as 
a coupling to a thermal reservoir in this paper. 

The paper is laid out as follows. The physical system is described in Section |H] along with 
the master equation describing the coupling to the reservoirs. The error correction protocol 
with measurement is analyzed in Sec. IIHI In Sec. IIVI we describe error correction without 
measurement and make a comparison with the continuous error correction described in 
followed by a discussion and conclusions in Sec. |V| 
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II. THE PHYSICAL SYSTEM 



We model a physical implementation of the three-bit repetition code protecting against 
bit-flip errors consisting of six individually addressable qubits, perhaps contained in an ion 
trap or array of ion microtraps. For simplicity we assume that the logical |0) and |1) states 
are degenerate and coherent transitions are achieved by two-photon Raman pulses. 

The accumulation of bit-flip errors in the qubits is modeled by a coupling to a hot thermal 
reservoir. In practice errors arise due to a noisy environment and imperfect control fields 
applied to the qubits. Here we simply replace these specific sources of errors with a coupling 
to a thermal reservoir with an effective temperature. The important parameter is the heating 
rate 7^ at which each qubit accumulates errors. For qubits with degenerate logical states 
|0) and |1) such a coupling gives rise to a completely mixed steady state. Quantum error 
correction slows the rate at which errors accumulate. 

The interaction Hamiltonian for the system is 

Hi = ( fi ( f ) a ™ + 9i (t)°i* + h i (*) Hi + Y^ % (*) P a } (!) 

8=1 I j>i ) 

where Pauli pseudo-spin operator acting on the zth qubit, Hi is the Hadamard 

operator acting on the ith qubit, and P^- is the two-qubit phase shift operator 

Pij = diag{a 00 , a 01 , a w , a n } (2) 

represented in the computational subspace of qubits % and j. This type of operator can 
be realized for example by applying a state dependent force on a pair of ions in separate 



microtraps 
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22|. We make the simplifying assumption that the pushing gate can be 
carried out between any pair of qubits without extra swap gates. This requires the ability 
to physically shuffle the qubits into proximity as needed. The parameters fi{t),gi(t),hi(t), 
and ki(t) represent applied fields which can be turned on and off at will. 

In addition to the unitary evolution the qubits are also coupled to two external thermal 
reservoirs at temperatures Th and T c . The reservoir at Th heats the qubits causing bit flips 
at a rate 7^ while the reservoir at T c cools the ancilla qubits at a rate T c . We assume that 
the cooling can be turned on and off at will either by lifting the degeneracy of the lo gica l 
states via an external field or by an active cooling method similar to that described in [18]. 
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FIG. 1: Quantum circuit for one round of error correction. 
The system evolution is determined by the master equation 
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where p(t) is an additional control parameter which allows a variable coupling to the cold 
reservoir and n c = (exp(^,a;//csT c ) — l) - is the thermal photon number of the cold reservoir 
at the energy separating the |0) and |1) states of the qubits. Note that this assumes the 
qubits are at least temporarily made nondegenerate. For an active cooling scenario the 
parameter n c characterizes the residual probability to find an ancilla qubit in the |1) state 
following the cooling step. 



III. ERROR CORRECTION WITH MEASUREMENT 

The basic quantum error correction protocol is illustrated in Fig. ^ It consists of a 
cooling step, ancilla preparation, transversal CNOT between the data qubits and ancilla 
qubits, ancilla decoding, syndrome measurement, and appropriate correction. When this is 
successful it corrects bit-flip errors of weight 1 on the data qubits. 

The fundamental gates available in the Hamiltonian do not include the CNOT gate. 
Therefore we express the CNOT as a sequence of fundamental gates as shown in Fig. |21 
This leads to the full error correction circuit in Fig. El 

The error correction protocol in Fig.^is simulated by the method of quantum trajectories 
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231 ] . In this method the master equation is unraveled as a sum over pure states 

fl£c(*))(^REcCOI (4) 

REC 

conditioned on a measurement record REC. The system undergoes a continuous evolution 
governed by a non-Hermitian Hamiltonian 

f r 6 r 6 \ 

H cS (t) = Hj(t) -i\Y + P(*)y K + !) S + P(*)y n " 2 (5) 

\ ~ i=4 i=4 / 

with discontinous bit-flip errors generated by the action of Oi X occurring with probability jhdt 
on all the qubits, heating events generated by <7; + occurring with probability p{t)T c n c (ai-ai + ) 
and cooling events generated by ov with probability p{t)T c {n c + l)(o" i+ o' i _) acting on the 
ancilla qubits. 

For simplicity we assume that the protocol takes place in 16 steps each of duration r. 
The gates are carried out in sequence by setting the fields fi(t), gi(t), hi(t), and fcj(i) to 
appropriate values, typically tt/2t. Note that the performance of the error correcting code 
depends directly on r since the probability for a bit-flip error is 7^r per qubit. 

A. Numerical results 

First we consider the data fidelity following each round of error correction. This is plotted 
in Fig. Ufa) and (b) for \ip{0)) = \0l)- in all cases the fidelity decays in time with the rate 
and steady state value depending on the quality of the ancilla preparation. The heating 
rate is 10~ 3 for all curves in the plot. One might at first expect the steady state value of 
the fidelity to be 1/8 for the 3-bit code because there are 2 3 = 8 states in the Hilbert space 
of the physical data qubits. This assumes that the data qubits will explore all of the basis 
states with equal probability. However, the error correction protocol drives the data qubits 



91 



Z(0 C ) 



P 



\Q2) ~ Z, 



X_ - Z(0 t ) -p-x + -z_ 



FIG. 2: Quantum circuit for a CNOT gate. Here 9 C = (aoo — ^lo) an d Ot = («oo — «oi)- The 
parameters are the diagonal elements of the pushing gate P. Z{9) is a rotation about the z-axis 
by an angle 6 and likewise for X{9). Z± = Z(±n/2) and likewise for X±. 
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FIG. 3: Quantum circuit for one round of error correction expressed in terms of operation available 
in the Hamiltonian. The gates and parameters are as defined in Fig. El 

preferentially to the states 1 000) and 1 1 11) even if many errors have accumulated in the data 
qubits. For very strong error correction the data qubits will effectively switch between 1 000) 
and | 111) giving F 2 = 1/2 while for weak error correction the data qubits will explore all the 
basis states equally giving F 2 = 1/8. Here strong error correction means that the protocol 
is likely to proceed without error for any single round while weak error correction means it 
is unlikely to proceed without error. 

The effect of the heating rate on the data fidelity is plotted in Fig. El^c). We have chosen 
values for 7 C and n c which give reasonable ancilla cooling. As the heating rate is increased 
the errors accumulate in the data qubits more rapidly while at the same time the ancilla 
qubits are less likely to be properly prepared before they are coupled to the data qubits.. 
For very strong heating, 7^ = 0.1, the fidelity is essentially constant at 1/8. 

We also consider the ancillla fidelity immediately following the cooling step. Several 
examples are plotted in Fig. El^d). These curves split into two groups. In the first (i, iii, 
iv) the ancilla fidelity is independent of the number of rounds preceding. This is the case 
whenever the cooling rate T c is large enough (> 3) that effective cooling takes place. Then 
the value of the ancilla fidelity is determined by n c or alternatively the temperature T c . The 
second group of curves (ii and v) show a decay of the fidelity in time to a steady state value. 
This occurs when the cooling rate is too small. In this case the cooling is not effective and 
the fidelity decays as errors accumulate in the ancilla qubits, either directly or by coupling 
to the data qubits. 
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FIG. 4: (a) Data fidelity for lh = 10" 3 , T c = 3 and n c = (i), 1(T 3 (ii), 1(T 2 (iii), KT 1 (iv). (b) 
Data fidelity for n c = 1(T 2 , j h = 1(T 3 , and T c = 30 (i), 3 (ii), 1 (iii), KT 1 (iv). (c) Data fidelity 
for n c = 10~ 2 , T c = 3, and = 10~ 3 x 1, 2, 3, 5, 10, 100 from top curve to bottom curve, (d) 
Ancilla fidelity for j h = 10" 3 and [n c ,T c ] = [10~ 2 ,3] (i), [10~ 2 ,1] (ii), [10~\3] (iii), [0.5,3] (iv), 

[io-Vo- 1 ] (v). 



B. Ancilla fidelity 



We have been able to understand some of these behaviors with simple models. For the 
ancilla cooling there are two types of extreme behavior to consider — fast (lines i, iii, and iv 
in Fig. HJl) and slow (line v in Fig. H]i) cooling. For fast cooling we assume the ancilla always 
reaches it's steady state during the cooling step. The populations under cooling follow a set 
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of rate equations 



P = AP l + AP 2 + AP 4 - 3BP 
Pi = AP 3 + AP 5 - (A + 2P)Pi + PP 
P 2 = AP 3 + AP 6 - (A + 2P)P 2 + BP 
P 3 = AP 7 - (2A + P)P 3 + BP 1 + BP 2 
P 4 = AP 5 + AP 6 - (A + 2P)P 4 + PP 
P 5 = AP 7 - (2A + P)P 5 + BP 4 + PPx 
P 6 = AP 7 - (2A + P)P 6 + BP A + PP 2 
P 7 = -3AP 7 + PP 6 + PP 5 + PP 3 



where Pj is the population of the ith ancilla state expressed as a decimal, A = T c ( 
and B = T c n c . In the steady state there is a symmetry P x = P 2 = P 4 = P a and P 3 
P 6 = P b which reduces these equations to the four 



= 3AP a -3BP 
= 2AP b -(A + 2B)P a + BP 
= AP 7 - (2A + B)P h + 2PP a 
= -3AP 7 + 3PP fe 



with the normalization condition 



P + 3P a + 3P 6 + P 7 = 1. 



The steady state solution is 




-3 



and the ancilla fidelity is 
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Note that the deviation from perfect fidelity is due solely to the nonzero temperature of the 
cold reservoir. Comparison with the simulation results shows good agreement. 

Now consider the case of slow cooling. If there is no cooling at all then the ancilla fidelity 
will decay because of heating of the ancilla qubits and the transfer of errors from the data 
qubits. At each measurement step the ancilla qubits are projected onto a single basis state. 
When there is effectively no cooling the fidelity is simply the probability that the ancilla is 
projected onto the |000) state. This is determined in turn by the fidelity of the data qubits. 
A comparison of the ancilla and data fidelities shows a close correlation for short times. For 
long times the ancilla qubits reach a different steady state than the data qubits. 

The steady state value of the ancilla qubits can be determined by a self-consistency 
condition from the dynamical equations governing the cooling. These equations are the same 
as given above. For simplicity we take B = and solve the equations for the individual 
initial conditions Pj = 1, P^ = since those are appropriate following the measurement of 
the ancilla state. The solutions are % — 0: 



Po(t) = 1, (11) 



i = l,2,4: 



Pi(t) = e- M (12a) 
P (t) = 1 - e~ At , (12b) 



3,5,6: 



Pi(t) = e~ 2At (13a) 

P hk {t) = e At - e' 2At (13b) 

P (t) = 1 - 2e' At + e~ 2At , (13c) 

where {i,j,k}= {3, 2, 1}, {5, 4, 1}, {6, 4, 2}, and i = 7: 

P 7 (t) = e~ 3At (14a) 

PsfiAt) = e- 2At -e- 3At (14b) 

Pi,2,4(f) = e- At - 2e- 2At + e- 3At (14c) 

P (t) = 1 - 3- At + 3e- 2At - e- 3At . (14d) 
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Defining x = e At we can write the fidelity as 

F = p +pi(l -x) +p 2 (l -xf +p s (l -xf (15) 

where pi is the probability to have % ancilla bits in state |1) prior to cooling. At steady 
state these are related to the fidelity. Roughly we can take p\ = (1 — F ss ) + F ss a where 
a — 16 steps x ^ h x 6 qubits = 967^ is the probability for the ancilla qubits to acquire an 
error during the round. Note that errors on the data qubits are mapped onto the ancilla 
qubits and contribute to the ancilla infidelity in the limit of no cooling. Putting F ss on both 
sides yields 

F ss = F ss (l -a) + [(1 - F ss )(l -a) + F ss a] (1 - x) (16) 
and solving for F ss yields 

F ss = ^-^-f . (17) 
1 — a — x + lax 

Again, comparison with the simulation results shows good agreement. 



C. Data fidelity 

The data fidelity following error correction can be found with a similar set of rate equa- 
tions, however the behavior is more complicated than the cooling of the ancilla. We trace 
the evolution of a set of probabilities Pi from one round of error correction to the next. 
There are 16 different events which can occur in one round of error correction. The ancilla 
can be properly cooled (1) or not (2), the ancilla can be properly prepared (1) or not (2), 
and we can have (1), 1 (2), 2 (3), or 3 (4) errors in the data qubits. We label these events 
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Ill, 112, . . ., 224 accordingly. The probabilities for these events to occur are 



Pill 


= F a [(l-af + a 3 ] 


(18a) 


Pll2 


= F a [(l-a) 3 + a 3 ] 3/?(l -/3) 2 


(18b) 


P113 


= F a [(l-af + a 3 } 3/3 2 (l - /?) 


(18c) 


PlU 


= F a [(1 - a) 3 + a 3 ] (3 3 


(18d) 


P121 


= F a [3a(l - a) 2 + 3a 2 (l - a)] (1 - /3) 3 


(18e) 


Pl22 


= F a [3a(l - a) 2 + 3a 2 (l - a)] 3(3(1 - (3) 2 


(18f) 


Pl23 


= F a [3a(l - a) 2 + 3a 2 (l - a)] 3/3 2 (l - (3) 


(18g) 


Pl24 


= F a [3a(l - a) 2 + 3« 2 (1 - a)] /5 3 


(18h) 


P211 


= P221 = (l-Fa)(l-/3) 3 


(18i) 


P212 


= P222 = (l-Fa)3/3(l-/3) 2 


(18j) 


P213 


= P223 = (1 - Fa)3/3 2 (1 - /?) 


(18k) 


P214 


= P224 = (1 - F a )/3 3 


(181) 



where F a is the ancilla fidelity following cooling, a = 157^ is the probability per qubit for 
a single error in the ancilla qubits during a round of error correction, and (3 = 167^ is the 
probability per qubit for an error in the data qubits during a round of error correction. Now 
we consider the flow of probability for each event. Given that there are data qubits in 
state |1) we have 

foo = Pin +P112 + 7^122 + 7^212 (19a) 
2 2 

foa = Pl21 + 7jPl23 +P211 + gP213 (19b) 

2 2 

fob = Pll3 + gPl22 + P124 + -P212 + P2U (19c) 
f07 = PlU + 7jPl23 + ^213- (19d) 

Given 1 data qubit in state |1) we have 

112 12 . . 

JaO = Pill + gPll2 + gPl21 + gPl23 + 77^211 + 7jP213 (20a) 

2 7 2 7 2 

faa = -P113 + gPl22 + ~Pl24 + ~P212 + ~P2U (20b) 

2 2 7 2 7 

fab = gPll2 + PlU + gPl21 + gPl23 + 7^211 + g^213 (20c) 

12 12 1 

fat = gPll3 + gPl22 + 7JP124 + gP212 + 7^214- (20d) 
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Given 2 or 3 data qubits we have the same equations as for 1 or respectively with Pq — > P-?, 
P a -> P b , P b -> P , and P 7 -> P . That is 



/&0 — /a7, 


(21a) 


/feci ,/a6 j 


(21b) 


J 66 J eta 5 


(21c) 


/&7 = /oOj 


(21d) 


/70 = fo7, 


(21e) 


fla = fob, 


(21f) 


fib = foa, 


(21g) 


fn — foo- 


(21h) 



This comes from the fact that the error correcting code treats 1 000) as equivalent to 1 1 1 1) . 

These have been presented as conditional probabilities, but they allow us to calculate the 
flow of probability from round i of error correction to round i + 1 when they are coupled 
with the actual probabilities from round i. This takes the form of a set of rate equations 



P (i + 1) = P (i) foo + P a (i)fao + Pb(t)fbo + Pi{i)ho, (22a) 

P a {l + 1) = Po{t)foa + Pa{i)faa + P b (l) fba + Pj^fla, (22b) 

P b {i + 1) = P (i)fob + Pa{i) fab + Pb(i)fbb + Fr(i)frb, (22c) 

P 7 {i + 1) = P (i)/o7 + Pad) far + P b (i)f b7 + P^frt, (22d) 



which we have solved numerically to produce plots in Fig.EJof the data fidelity for comparison 
with the simulations. 

It turns out that for the case of good ancilla cooling the decay of the fidelity is very nearly 
exponential following a relatively large decrease in the initial rounds of error correction. We 
consider these two features in turn. The large initial decay in the data fidelity can be found 
by considering Pq{1) from equation ()22a|) . 

Po(l) = foo=(^^J(l-3a-(3) + f3 (23) 

where we have used the definitions from Eqs. (10), (18), and (19) keeping only first order 
terms in a and (3. The form of the fidelity following one round of error correction reveals that 
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FIG. 5: Comparison of trajectory results and rate equations for n c = 10~ 2 , T c = 3, and 7^ = 0.001 
(i), 0.002 (ii), 0.1 (iii). 

the source of the large initial decay in the fidelity is a combination of non-zero temperature 
of the cold reservoir (n c 7^ 0), errors on the ancilla qubits, and errors on the data qubits. 
Either imperfect ancilla cooling or errors on the ancilla qubits during preparation will always 
contribute to a reduced data fidelity following the first round of error correction while errors 
in the data qubits only contribute in combination with imperfect ancilla cooling. 

We quantify the subsequent slow decay of the fidelity by solving the rate equations to a 
few orders in the heating rate a and matching term by term with decaying exponentials of 
the form 

Po(n) = P ,ss + (Po(k) - P , ss )So (n ~ k) (24a) 
Pain) = P a , ss - (P a , ss - P a (k))5~^ (24b) 
P b (n) = P b , ss - (P b , ss - P b (k))5- {n - k) (24c) 
Pr(n) = Pr,ss ~ (P 7 ,ss - P 7 (k))5^ n - k) (24d) 

and find the form of 5 to the same order in a where k is a small number to account for the 
large initial decay of the fidelity. Here we take f3 = a and F a — 1 for simplicity. To illustrate 
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this, consider the fidelity written to second order in a 



Po(n) = / n o + (n - l)/oo" 2 /oa/ao + 
= 1 - 3a + (33 - 21n)a 2 



(n-2)(n-3) 



Od 



(25) 



with the steady state value 



P< 



ho + f. 
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0,ss 



We write Sr, as 



2 /oi + /02 + fio + /l3 
- (l - 3a + 24a 2 ) . 



= ao (l + aia + 02a 2 ) 



(26) 



(27) 



where dj are coefficients to be determined. We choose k = 4 in order to keep all the terms 
in Po(n) to second order and write out Eq. I24al as 



P (n) = i - |a + 12a 2 + 



1 / 3\ , 2 v 2 1 

- + I 2ai I a + (2a 2 + 3a x — 6a± — 63) a a\ 

2 V 2 y 2 



raa 



3 7q 1\ o Iqqo 
+ ( -a a - -a x - -a 2 I na + -a^ a 



(28) 



Matching the constant term requires ao = 1, matching the a term requires ai = 0, and 
matching the na term requires a 2 = 42. Then we have 5 = 1 + 42a 2 . 



D. Cyclic features in error correction 

So far we have considered the gradual decay of the fidelity over many rounds of error 
correction. Now we look at the fidelity and the entropy at each step of the error correction 
protocol (see Fig. EJ). We use the von Neumann entropy defined as 

S = -tr(p]ogp) (29) 

where the density operator p is found by averaging It/'recWX'^recWI over many trajectories. 
These are plotted in Fig. El showing the cyclic process when the error correction is working, 
particularly for the ancilla qubits. We have chosen 7^ = 10~ 3 and two sets of parameters: 
{n c , 7c} — {10~ 2 , 3} (i) for effective error correction and {10 -2 , 10 -2 } (ii) for ineffective error 
correction. The difference is apparent in both the fidelity and the entropy. When the error 
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correction is effective there are periodic increases in the fidelity and corresponding decreases 
in the entropy which appear in the time step in which the correction is applied. When the 
error correction is ineffective there is simply a steady decrease in the fidelity and increase in 
the entropy of the data qubits. 

It is worthwhile considering the physical reasons for the differences between the sets 
of curves labeled (i) and (ii) in Fig. |U1 In the curves (i) the ancilla is well-cooled; after 
every error-correction cycle, its entropy goes down to the same "floor" value and the ancilla 
fidelity rises to the same "ceiling" . That this floor is not zero is due to the nonzero reservoir 
temperature represented by n c and the finite cooling rate T c . However, even though the 
ancilla entropy does not increase on average in this case, the data entropy continues to 
increase on average with small decreases upon successful error correction; this is due to the 
fact that there is always a probability that an uncorrectable error will take place (an effect 
proportional to 7^) . The nonzero ancilla entropy (and thus n c ) also plays a role in the steady 
increase of the data entropy because an incorrectly prepared ancilla will cause the wrong 
syndrome to be extracted so the data gets miscorrected. 

For the curves labeled (ii) there is effectively no cooling of the ancilla qubits and the 
ancilla entropy steadily increases. Without a properly prepared ancilla the syndrome is not 
extracted properly and the errors in the data qubits cannot be corrected. 



IV. ERROR CORRECTION WITHOUT MEASUREMENT 

It is also possible to implement error correction without measurement by means of ad- 
ditional quantum gates to carry out the appropriate correction conditionally based on the 
state of the ancilla qubits. A circuit for implementing the three bit repetition code using 
two ancilla qubits is shown in Fig. The correction is carried out by means of three Toffoli 
gates which flip the approriate data qubit based on the state of the ancilla qubits. 

We use the same Hamiltonian to simulate this quantum circuit with an appropriate 
sequence of control fields. In order to do this we must expand each Toffoli gate into a 



sequence of one- and two-bit gates. Based on results in Ref. [24J one such circuit is shown 
in Fig. |H1 The CNOT gates appearing in the error correction circuit and the Toffoli gates 
are again expanded according to Fig. El and a sequence of fundamental gates is found to 
implement the error correction without measurement in 68 equal time steps. 
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FIG. 6: Data fidelity (a), ancilla fidelity (b), total entropy (c), data entropy (d), and ancilla 
entropy (e) for [n c ,T c ,j h ] = [W~ 2 , 3, 10 -3 ] (i) and [1(T 2 , 10 -2 , 10 -3 ] (ii). In plot (b) the thin line 
corresponds to parameters (i) and the thick line corresponds to parameters (ii). 

It should be noted that this scheme is quite similar to that described in Ref. [19| but with 
three important differences. First, Sarovar and Milburn consider a continuous coupling of 
the ancilla qubits to a zero temperature reservoir while we assume that the ancilla qubits 
are cooled for a single time step at the beginning of the protocol by coupling to a finite 
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FIG. 7: The quantum circuit diagram for implementing the three qubit repetition code without 
measurement. 
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FIG. 8: The quantum circuit to implement a Toffoli gate using one- and two-bit gates. 

temperature reservoir. Second, we include heating of the ancilla qubits as well as the data 
qubits. Third, they use a time independent Hamiltonian to carry out the error correction 
but it is not clear how that might be realized experimentally. We have used an explicit 
sequence of gates with a time dependent Hamiltonian which could be implemented in a 
variety of physical systems. 

First we consider the data fidelity following each round of error correction plotted in Fig. El 
for |^(0)) = \0l). In all cases the fidelity decays in time with the rate and steady state value 
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FIG. 9: (a) Data fidelity for j h = 10" 4 , T c = 3 and n c = (i), 1(T 3 (ii), 1(T 2 (iii), KT 1 (iv). (b) 
Data fidelity for n c = 10~ 2 , T c = 3, and j h = 1(T 4 x 1, 2, 3, 5, 10, 100 from top curve to bottom 
curve. 
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depending on the quality of the ancilla preparation. The heating rate is 10~ 4 for all curves 
in the plot (a). Just as for error correction with measurement, one might expect the steady 
state value of the fidelity to be 1/8 for the 3-bit code because there are 2 3 = 8 states in the 
Hilbert space of the physical data qubits. However, for very strong error correction the data 
qubits will effectively switch between 1 000) and |111) giving F 2 = 1/2 while for weak error 
correction the data qubits will explore all the basis states equally giving F 2 = 1/8. 

The effect of the heating rate on the data fidelity is plotted in Fig. HJb). We have chosen 
values for 7 C and n c which give reasonable ancilla cooling. As the heating rate is increased 
the errors accumulate in the data qubits more rapidly while at the same time the ancilla 
qubits are less likely to be properly prepared. For very strong heating, 7^ = 0.01, the fidelity 
is essentially constant at 1/8. Notice that 7^ is an order of magnitude smaller in these plots 
that for error correction with measurement shown in Fig. 0] This is necessary because it 
takes many time steps to carry out the Toffoli gates giving a larger failure rate per round of 
error correction. 

The fidelity and entropy for error correction without measurement is shown in Fig.llOlfor 
each time step. Here the ancilla qubits have a strong cyclic behavior while the data qubits 
show an overall increase in entropy along with cyclic variations. 

V. DISCUSSION AND CONCLUSIONS 

Our analysis has revealed a similar cyclic pumping of entropy from the data qubits to the 
ancilla qubits to the cold reservoir in for error correction with and without measurement. 
They also share a similar decay of the fidelity of the data qubits due to incomplete removal 
of the entropy with a fast initial decay followed by a slower decay as the total system reaches 
a steady state. The performance is quite sensitive to the ancilla cooling with either slow 
cooling or cooling to a nonzero temperature leading to degraded performance of the error 
correction protocol. It is interesting to note that in these circumstances the decay of the 
fidelity is first order in the heating rate as a single bit flip in the ancilla qubits is enough to 
disrupt the error correction protocol. 

In fault-tolerant quantum error correction schemes (see Ref. [5] for example) the linear 
dependence on the error rate is addressed in two ways. The ancilla state is verified prior 
to coupling to the data qubits leading to a reduction in the probability that an incorrectly 
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FIG. 10: Data fidelity (a), ancilla fidelity (b), total entropy (c), data entropy (d), and ancilla 
entropy (e) for [n c ,T c ,j h ] = [1(T 2 , 3, 10" 4 ]. 

prepared ancilla state leads to an incorrect syndrome. Furthermore, multiple syndrome ex- 
traction leads to a reduction in the probability that the error in the data qubits is incorrectly 
diagnosed and therefore miscorrected. 

It should also be noted that the ancilla cooling imposes a minimum energy requirement 
(this is separate from the minimum energy requirements for quantum logic discussed by 



one of us in 



17| ) to perform error correction due to Landauer's priniciple which requires 
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/cbT c In 2 of energy per bit of information erased. 
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